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Abstract: The state of spatially correlated down-converted photons is 
usually treated as a two-mode Gaussian entangled state. While intuitively 
this seems to be reasonable, it is known that new structures in the spatial 
distributions of these photons can be observed when the phase-matching 
conditions are properly taken into account. Here, we study how the vari- 
ances of the near- and far-field conditional probabilities are affected by the 
phase-matching functions, and we analyze the role of the EPR-criterion 
regarding the non-Gaussianity and entanglement detection of the spatial 
two-photon state of spontaneous parametric down-conversion (SPDC). 
Then we introduce a statistical measure, based on the negentmpy of the joint 
distributions at the near- and far-field planes, which allows for the quantifi- 
cation of the non-Gaussianity of this state. This measure of non-Gaussianity 
requires only the measurement of the diagonal covariance sub-matrices, and 
will be relevant for new applications of the spatial correlation of SPDC in 
CV quantum information processing. 
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1. Introduction 

In spontaneous parametric down-conversion (SPDC), photon pairs are generated with several 
degrees of freedom quantum correlated. In particular, the down-converted photons are spatially 
entangled. Due to energy and momentum conservation, the sum of the transverse momenta and 
the difference of the transverse positions of the photons can be well defined, even though the 
position and momentum of each photon are undefined [1,2], This type of Einstein-Podolsky- 



Rosen (EPR) correlation [3], has been used as a resource for fundamental studies of quantum 
mechanics [1,2,4-6], for quantum imaging [7, 8], and experiments of quantum information 



The process of SPDC has been studied extensively in the past [ 12-14], and much of the recent 
effort has been put in the quantification of the spatial entanglement for a given experimental 
geometry. Traditionally, this has been done through the technique of Schmidt decomposition 
of the two-photon wave function, which gives the Schmidt number, K, and the Schmidt modes 
allowed for each photon [15-18]. While this approach describes important properties of the 
spatial correlation, it basically gives no information of the form of the spatial joint distributions, 
therefore giving no information about the Gaussianity of the spatial two-photon state of SPDC. 

Moreover, it is well known that the state of the down-converted photons depends on the 
phase-matching conditions. Nevertheless, due to its complex structure, the phase-matching 
functions are usually approximated by Gaussian functions [17-22]. While the approach be- 
ing adopted seems to be reasonable, it has already been shown that fine (new) structures in the 
spatial distributions of these photons can be observed due to (the manipulation of) the phase- 
matching conditions [23]. 

In this work we study the non-Gaussianity of the spatial two-photon state of SPDC by prop- 
erly taking into account the phase-matching conditions. We start by showing how the vari- 
ances of the near- and far-field conditional probability distributions are affected by the phase- 
matching functions. Then, we analyze the role of the EPR-criterion [3, 24] regarding the non- 
Gaussianity and entanglement detection of this state. Even though it has been demonstrated 
that higher order separability criteria can be used for the entanglement detection of spatial non- 
Gaussian entangled states [25], we show that a proper consideration of the phase-matching 
function reveals, precisely, when the simpler EPR-criterion can still be used for the spatial en- 
tanglement detection. We also show that (and when) the EPR-criterion can be used as a witness 
for the non-Gaussianity of this state. 

Furthermore, we introduce a statistical measure, based on the negentwpy [26] of the near- 
and far-field joint distributions, which allows for the quantification of the non-Gaussianity of the 
spatial two-photon state of SPDC. This measure does not correspond to a quantum mechanical 
generalization of the negentropy, such as the non-Gaussianity measure based on the quantum 
relative entropy (QRE) [27, 28], and so does not require the knowledge of the full density 
matrix. Only the moments associated with the diagonal sub-matrices of the covariance matrix 
need to be measured. Thus, it is experimentally more accessible [29]. Moreover, for most of the 
configurations used so far, we show that its value can be estimated from the (easier to measure) 
marginal and conditional distributions. We also demonstrate that it has common properties 
with previous introduced measures of non-Gaussianity [27,28]. The quantification of the non- 
Gaussianity of a quantum state has important applications for quantum information [25, 27, 
28], and thus the practicality our measure shall be relevant for new appUcations of the spatial 
correlations of SPDC in this field. 

2. The phase-matching conditions and the variances of the conditional probabilities 

We consider the process of quasi-monochromatic SPDC, in the paraxial regime, for configu- 
rations with negligible Poynting vector transverse walk-off, that can be obtained using non- 
critical phase-matching techniques [23]. In the momentum representation, the two-photon state 
is given by [12-14] 



[9-11]. 




(1) 



where 1 9i , ^2 ) represents a two-photon state in plane-wave modes whose transverse wave vec- 
tors are qi and qT- Lis the crystal length, kp is the pump beam wave number, ^ ~ Wq~ 2iz/kp, 
and Wo is the pump beam waist, which is located at z = 0. This state may be rewritten in the 
coordinate space as [23, 30] 



(2) 



where we define the function sint as sint (x) = ^ dt sine f = 1 — -I Si (x). Si [x] being the sine 
integral function. The functions sine {b\q\^ /2) and ^ sint ( \p\'^ /2b) form a Fourier transform 
pair 

The sine and sint functions arise from the phase-matching conditions, and due to the dif- 
ficulty of deaUng analytically with them, they are usually approximated by Gaussian func- 
tions. The approach that has been adopted consists in approximating the function sine [bx^) by 
^-abx^ ^ Sometimes it is used that a = 1 [15, 19], and in other cases the value of a is chosen 
such that both functions coincide at [17, 20, 21] or at 1/e [22] of their peak. While this 
approximation seems to be reasonable for the entanglement quantification [15], there has been 
no investigation to determine how precise it is for describing the distribution of the momentum 
correlations. Besides, it should be noticed that once a Gaussian approximation is adopted for 
the sine function, the corresponding approximation for the sint function is already defined by 
the Fourier transform. Therefore, it is also not clear that such approximation is indeed good for 
describing the position distributions. Thus, it is not clear whether the SPDC two-photon state 
can indeed be written as a two-mode Gaussian state. 

To investigate this point we study how the variances of the momentum (far-field) and position 
(near-field) conditional probability distributions are affected by the phase-matching function, 
and compare the obtained results with the cases where Gaussian approximations are considered. 
Due to the symmetry of the two-photon wave functions, there is no loss of generalization if we 
work in one dimension (i.e., y\ —yo—Q and qy\ = (7,2 = 0). To simplify our analysis we define 
the following dimensionless variables: xj = xj/wo and qj ~ woqj, j ~ 1,2. The probabilities 
for coincidence detection at the far- and near-field planes are 



p'pAquqi) - I sine - ---- ) , (3) 



pU^U^2) - .-5^(-+-'^(sintl£i^)', (4) 



where = {zq +z^)/zq, P = yjLj (2zo), and zo = ^pW'o/2 is the diffraction length of the 
pump beam. The joint probabilities are related with the conditional (and marginal) probabil- 
ities p%p{q\ \q2) \.p%F{qi)\ plip{xx \x2) [p^f (iy)], through the rule: p{^y\^2) = 

In the case of the Gaussian approximations discussed above, the probabilities of coincidence 
detection are 

P^.,„,(ii,i2) - e-^^^^^^^e-^''^-'^^, (6) 

where different values of a, represent distinct Gaussian approximations for the sine function. In 
Figure 1(a) [(b)] we compare the curves p^p and Ppp^. (p^p and p'^pa ) considering i2, ^2 = 0, 
a = 1 (and the crystal centered at z = 0), for the case where the dimensionless parameter 




Fig. 1. The far- and near-field conditional probability distributions, while considering the 
state of SPDC (Pp-f and p^ip) and when some Gaussian approximations are assumed 
{ppfa and p^fg, )- In {a)-{d) the curves p^pp and pf^p (red solid lines) are compared with 

p'ff a2 ^"'^ PNFai ^"-"^ '^i ~ ^'^^ (green dashed line), a2 = 0.72 (blue dot-dashed line), and 
for some specific values of P. In (e) and (/) the variances of the momentum and position 
conditional probabilities are plotted in terms of P, respectively. a3 = 1 . 



P = 0. 1 . This parameter has been used in the study of the quantification of the spatial entangle- 
ment [15, 16], and it brings universality to the theory since a certain value of P can be reached 
in three different ways. Here we considered the values of L and kp as fixed parameters such that 
P varies with wq. In this case we find ai = 0.45 (a2 = 0.72) for the case where the sine and 
Gaussian functions coincide at 1/e (1/e'). From Figure l(a)-(b) one can see that for a small 
value of P, the Gaussian approximation only describes properly the momentum conditional 
distribution. The position conditional distribution is barely described by the approximation. In 
Figure 1(c) and Figure 1(d) we have the same type of analysis but now for a larger value of 
P. In this case, the Gaussian approximation is useful only for describing the position condi- 
tional distribution. The overall behavior of the Gaussian approximations is showed in Figure 
1(e) and Figure 1(f), where the normalized variances of the far- [{Aqi \q^)^L/kp] and near-field 
[{Axi \x2f'kp/L\ conditional distributions are plotted in terms of P. 



3. The EPR-Criterion as a witness for the non-Gaussianity of the spatial two-photon 
state 



The near and far-field conditional probabilities can be used for implementing the EPR-paradox 
[1, 3,21, 24]. This is done by observing the violation of the inequality {Axi | „ li/^)^ > \, 

which certifies the quantum nature of the spatial correlations of the down-converted photons. 
Since we have determined how the phase-matching function affects the variances (Axilv,)^ 
and (Aq'il^^)^' we can also look for its effect on the EPR-criterion. This is showed with the 
red (solid) line in Figure 2, which was calculated for the values of X2 and q2 at the origin. 
For smaller values of P, the conditional variances are independent of the X2 and q2 values [1]. 
Whenever P increases, the variances become dependent on X2 and q2- Nevertheless for smaller 
values of X2 and q2, which are of most experimental relevance, the red (solid) curve shown in 
Figure 2 captures the overall behavior of the product of (Axi and (A^i for the state of 
Eq (1). As we can see, for values of P smaller than 0.56 or greater than 2.58, the EPR-criterion 
can safely be used for detecting the spatial entanglement of the two-photon state of SPDC. 

Besides of being useful as a entanglement witness, the EPR-criterion can also be used as 
a witness for the non-Gaussianity of the spatial state of SPDC. This emphasizes another ap- 
plication for this criterion, which has been related already with other quantum information 
tasks [24]. To observe this, we note that for a pure two-mode Gaussian state it is possible 
to show that (Axi \x2)^{Aqi = [see Appendix A], such that | is an upper bound for the 
EPR-criterion with these states. Since it has been demonstrated in [15] that the two-photon state 
of Eq. (1) is always entangled, we can say whenever the product of variances is greater than 
J, that it witnesses the non-Gaussianity of the entangled spatial two-photon state of SPDC. As 
one can see in Figure 2, this happens for 0.56 < P < 2.58. The Schmidt decomposition of the 
state wave function, used together with the EPR-criterion value, reveals the non-Gaussianity 
of a two-mode entangled state. Such observation does not necessarily hold true when other 
second-order moments criteria are considered. This is shown in Appendix B for the criterion of 
Ref. [31], and for the states considered in Figure 2. 
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Fig. 2. The EPR-criterion plotted in terms of P. The red (solid) curve corresponds to the 
EPR values of the two-photon state generated in the SPDC [Eq. (1)]. The green (dotted), 
blue (dot-dashed) and black (dotted) curves describe the EPR-criterion for the two-photon 
Gaussian states defined in terms of ai , a2 and 0:3 , respectively. 



4. Quantifying tlie non-Gaussianity of tlie spatial two-photon state of SPDC 

From our previous analysis it is clear that the spatial state of SPDC can not be seen as a two- 
mode Gaussian entangled state, even when it is generated with a Gaussian pump beam and in 
the case of perfect phase-matching. We now proceed to quantifying the non-Gaussianity of this 
state. First, we introduce the concept of negentropy which is the base of our approach [26]. 
The negentropy of a probabiUty density function p(^i,^2) is defined as N = H[p^{£,i,£,2)] — 
H[p{^i,S,2)], where p'^{^i,£,2) is a Gaussian distribution with the same expected values and 
covariance matrix of p{S,i , ^2)- The function H[p{^i , ^2)], called differential entropy, is defined 
asH[p{^i,^2)] = - .[d^id^2P{^i,^2)iog2p{^i,^2) [32]. The advantage of using negentropy is 
that it can be seen as the optimal estimator of non-Gaussianity, as far as density probabilities 
are involved. This is due to the properties that it is always non-negative, and that it is zero only 
for Gaussian distributions. Besides, it is invariant under invertible linear transformations [33]. 

Motivated by these properties we define the total non-Gaussianity of the spatial two-photon 
state of SPDC as 

nG^=N[p'pp{quq2)]+Nljy'^pixuX2)], (7) 

where N[ppp{qi,q2)] and N\pf;p{xi,X2)] are the negentropies of the far- and near-field joint 
distributions of Eq. (3) and Eq. (4), respectively. In Appendix C we give, explicitly, the cal- 
culations of N[ppp{qi,q2)] and N[pf^p{x[,X2)]- It is possible to observe that nG^ = if and 
only if \W) is a two-mode Gaussian state [see Appendix D]. Otherwise nO^ > 0. We obtain 
that N\ppp{qi,q2)] ~ 0.15 and that N[pfip{xi ,X2)] ~ 0.22. Thus, the total non-Gaussianity of 
the spatial two-photon state of SPDC [Eq (1)] is nG^ « 0.37. It is interesting to note that it 
does not depend on P. This was expected since the phase-matching functions do not change 
their functional form when P varies. The calculations performed for these negentropies can 
be adapted for different experimental geometries, or used for the proper determination of other 
quantum information quantities related with the differential entropy of the spatial joint distribu- 
tions [34]. The fact that nG^ 7^ 0, emphasizes that spatial Gaussian approximations should be 
taken carefully due to the extremality of Gaussian states [35]. For comparison purposes, we 
calculated in Appendix E the value of 5b, which is the measure of non-Gaussianity based on 
the QRE [27,28]. We obtain that 8b = 108 and such result also highlight the non-Gaussian 
character of the state of Eq. (1). Furthermore, it has the same behavior of nG^ , since it does not 
depend on the parameter P. 

In analogy with Eq. (7), we define the non-Gaussianity of the conditional and marginal dis- 
tributions as: nG^ = A?[pf-f fel^j)] +N[plip{xi\xj)] and nG^ = A?[/4f (*)] +N[plip{xi)] with 
ij = 1,2 and ; 7^ j. According to these definitions, one can observe that the non-Gaussianity 
of the spatial state of SPDC decreases under partial trace, such that nG^ > nG^; and that it 
is additive when the composite system is represented by a product state, i.e., if | ^) is a prod- 
uct state, then nG^ = 2nG^ [see Appendix D]. These are common properties with the QRE 
measure of [27,28]. 

In Figure 3(a) we plot the negentropies N[ppp {qi \q2)] and N[pfjp {x\ \x2)] as a function of P. 
It is interesting to note that these curves quantify the idea already presented in Figure l(e)-(f). 
As larger the value of P is, the less the conditional momentum distribution can be approximated 
by a Gaussian function. On the other hand, the conditional position probabilities tends to a 
normal distribution when P increases. In Figure 3(b) we plot the negentropies of the near- 
and far-field marginal probabilities. One can see that they have a different dependence on P 
in comparison with the conditional probabilities. Now, the near-field distribution tends to a 
Gaussian function for larger values of P, and the far-field one for smaller values of P. In Figure 
3(c) and Figure 3(d) we have nG'' and nG^ plotted in terms of P. The insets of these figures 
show the corresponding near- and far-field distributions at the points of minimum, indicated 
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Fig. 3. In (a)[(b)] the negentropies of the near- and far-field conditional (marginal) distribu- 
tions are plotted in terms of P. In (c) [(d)] we show the non-Gaussianity of the conditional 
(marginal) distributions. 



with red circles. 

As it is shown in Appendix F [See also Figure 3(c) and Figure 3(d)], in the limit of P ^ 1, 
one can use the relation between the joint and conditional density probabilities to decompose 
nG^ as the sum of iiG'^ and nG^: 

« + nC^ (iff P < 1 ) . (8) 

In a typical experimental configuration for SPDC, where the pump beam spot size is around 
1 mm at the crystal plane, the value of P can be smaller than 0.05. Thus, in general, the total 
non-Gaussianity of the spatial two-photon state of SPDC can be estimated in terms of the (easier 
to measure) near- and far-field negentropies of the conditional and marginal distributions. This 
simpUfy the measurement of nG^ , since there is no need to scan the whole transverse planes 
associated with the near- and far-field joint distributions. 

5. Conclusion 

We have investigated the spatial distributions of the entangled down-converted photons by 
proper considering the phase-matching conditions. By understanding how the near- and far- 
field conditional distributions are affected by the phase-matching function, we could show that 
the EPR-criterion [3] can be used as a witness for the non-Gaussianity of the spatial state of the 
SPDC. The work culminated in the quantification of the non-Gaussianity of this state, which 
was based in a new and very experimentally accessible measure. As it has been discussed 
in [27,28], the quantification of the non-Gaussianity of a quantum state has many applications 
in the area of continuous variables quantum information processing. Thus, we envisage the use 
of our result for new applications of the spatial correlations of SPDC. 



A. EPR-criterion for spatial Gaussian two-photon states 



As we mentioned in the main paper, the momentum representation of the spatial two-photon 
state generated in the spontaneous parametric down-conversion (SPDC) process, under paraxial 
approximation and for configurations with negligible Poynting vector transverse walk-off, can 
be written as [12-14] 

I'P) = fdqi f dq2'V{qi,q2)\lqi)\lq2), (9) 



where the amplitude , ^2) is given by 



(10) 



with A{qi + q2) being the angular spectrum of the pump beam and &{qi— q2) representing the 
phase-matching conditions of the non-linear process. If both functions are represented by Gaus- 



exp 



sian functions of the form h.{qi+q2 
[15, 17-22], we have in momentum representation that 

and in position representation that 



and 0(^1—^2) — exp 



2Si 



(11) 



(12) 



The parameter 5 can be adjusted in order to approximate the phase-matching function by 
distinct Gaussian functions. For a pump laser beam with a Gaussian transverse profile, = 
2/c^, where c is the radius of this pump at the plane of the non-linear crystal. 

Based on Eqs. (11) and (12), we can obtain the probability density functions for the condi- 
tional position and momentum distributions of the down-converted photons. The variances of 
these curves are, in general, independent of the value considered forx2 and ^2 [1] and here we 
use, for simplicity, X2 = and q2 = 0. In this case P{xi \x2 = 0) and P{qi \q2 = 0) are given by 



and 



P(xi|x2 = 0) 



Piqi\q2 = 0) 



51 



20r 



-exp 



err 



52 



— 9 



-exp 



51 



(13) 



(14) 



^al5l 

Here C7^ and 5l are the widths of the respective Gaussian functions [15]. The variances of the 
conditional distributions can be calculated directly from Eqs. (13) and (14), which results in 



and 



(^l|x2)" 



Then, the product of them is given by 

{AxiU,f{Aq,\,,) 



al + 51' 

al5^ 
2ial + 5l)' 

al5l 



{al + 5lf 



(15) 



(16) 



(17) 



Let us now consider the Schmidt number K for Gaussian states. It was showed in Ref. [15] 
that it is given by the expression 

Comparing Eqs. (17) and (18), it follows immediately that the product of the conditional vari- 
ances is a function of the Schmidt number, such that 

This expression is valid for any Gaussian approximation taken for the phase-matching func- 
tion. It is possible to see that the maximal value of {Axi | „ )^ (A^i )^ is equal to ^ and that it 
happens when the Schmidt number K = I (i.e., for product Gaussian states). 

B. Mancini et al. Criterion for the spatial entanglement of SPDC 

Another entanglement detection criterion, based on second-order moments, is the one intro- 
duced by Mancini et al. [31]. For the case of spatial entanglement it reads [29] 

[A(ii-i2)]'[A(^i+^2)]'> 1. (20) 

When this inequality is violated, the state is not separable. If this condition is satisfied, no 
information can be drawn. For the far and near-field joint spatial distributions [Eq. (3) and 
Eq. (4) of main paper, respectively], the variances in (20) are given by 

[A(^l +^2)]^ - j dq+qle-'^'i"^ = 1, (21) 

[A(ii ~X2)\^ °^ J dx^xl sint^^i^ ^ 4^P^ (22) 
where q+ ^ qi + qi and i_ = ii — X2- The constants 

Al = / sint'^2 ^ 1.4008, (23) 



and 



J<^<^^sint2,^2^ 0.5897. (24) 



Figure 4 shows the dependence of the Mancini-Criterion with the parameter P, while consid- 
ering the state of SPDC and when distinct Gaussian approximations for this state are considered. 
Since there is no upper limit for this criterion with Gaussian two-mode states, it is not possi- 
ble to use the Mancini-Criterion for the detection of the non-Gaussianity of the spatial state of 
SPDC. 

C. Negentropy of far- and near-field joint distributions of SPDC 

C.l. Negentropy of a probability density function 

Let p{^i,^2) be a probability density function. The Negentropy (A^) of p{^i,^2) is defined 
as [26] 

A^^//b^(^i,^2)]-//[p(^i,&)], (25) 

where p^{^i,^2) is a Gaussian distribution with the same expected value vector jU = {ji^^,!!^^} 
and same covariance matrix A of p{^i,^2)- The function , ^2)], called differential en- 

tropy, is defined as [32] 

H[pi^u^2)]^- fd^Jd^2Pi^u^2)iog2ipi^u^2))- (26) 



p 



Fig. 4. Comparison of thie Mancini-criterion while considering the state of SPDC (red sohd 
line) and when some Gaussian approximations are adopted. Here we use «! = 0.45 (green 
dashed line), ao = 0.72 (blue dot-dashed line), and = 1 (black dotted line) to describe 
the spatial Gaussian approximations [See Eq. (5) and Eq. (6) of main paper]. 



C.2. Negentropy of far-field joint distribution 

The joint probability density function in the far-field plane is given by 

p^pp{qi,q2) = ^-j^e 2 sine —{q\~q2) ■ (27) 
The eo variance matrix of Pff{qi,q2) is given by 

The elements in the diagonal are the variances of the marginal distributions Pff{qi) and 
Ppp{q2)- In order to calculate the terms of A^^, we first note that the expected values 
Mff = {(^i)i (^2)} are null. So, we have that 

(A^i)l = (A^2)| = JJ dqidq2q]pFF{quq2) 
{qiqi) S = {qiqx) S = JJ dqidq2qiq2PFFiqhq2) 

(30) 



+ (29) 



and 



4 V P^, 

Let us consider the following Gaussian approximation for describing the joint probability 
density distribution at the far-field plane: 

pUquq2)^^_e-^e-^. (31) 



The elements of the covariance matrix A^^ of ppp (^i , ^2) are given by 



{Aqif(y = {Aq2)c = j j dqidqiqlpFpiquqi) 



(32) 



and 



{qiqijG = {q2q\)G 



j j dqidqiqiqiPppiqi^qi) 



From Eq. (29), (30), (32) and (33) we can observe that the condition 

3 



5i = 



p2' 



(33) 



(34) 



implies that the joint distributions (27) and (31) have the same covariance matrix. Hereafter, 
the Gaussian distribution with the same covariance matrix of (27) is denoted by p'^fiqi^qi)- 
It is possible to show that the differential entropy of the Gaussian distribution Pff{qi,q2) is 



HlPFFiquq2)]^^og2 

The differential entropy of Pff{qi,q2) is given by 



ne 



V3 



(35) 



H[pFFiquq2)] 



log2 



J J dqidq2PFFiquq2) logi (Pff 1 , ^2 ) ) 
A-VlK 3P 



3P 4V2;r 



dq\dq2 j e 

p2 



Hi +5/2)' ji J 

logje 2 +log2smc —(^1-^2) 



(36) 



Making 



g M = and V = ^ (^1 — ^2) we obtain that 



^^[Pff (^1,^2)] = log2 



4V27r 1 
3P 21n2 

1 r= / c/y sine ^v^ In (sine ^v^"* 

v/TT Jo ^ 



(37) 



V sine ^ In ( sine ^ ) 



The last integral in (37) must be calculated numerically and it gives that 

-0.364. 

Therefore, the differential entropy of (^1,^2) is 

H[pfF{qhq2)] « iog2 + 1.17. 



(38) 



According to Eq. (25), we have that 

N[p^{q\A2)] = H[p'^p{quq2)]-H[l4p{qi,q2)] 

. 7teV3 \s/2% 
~log2— ^ log2^^-117 

/27 

wO.15. (39) 
C. 5. Negentropy of near-field joint distribution 

The joint probability density function in the near-field plane is given by 

p%p{xuX2) = -^e ^sint^^i^^, (40) 

where C is a normalization constant given by 

C=V2nAy^P. (41) 
The covariance matrix for the near-field joint distribution is 

s _( {x\X2)s 



where 



Kf = \ / n ;v \2 , (42) 



(Ail)! = (Ai2)5 = j j dx\dx2x\pljp{x\,X2) 



1 / 4A2P^ , 



and 



{XiX2) s ^ {X2Xi) s = 



J I dxxdx2XiX2pljp{x\,X2) 



Let us now consider the following Gaussian approximation for the joint probability distribu- 
tion at the near-field plane: 

P%(xx,X2)^^^e ^e . (45) 

The covariance matrix A^^ of p'^p ixi^xj) is determined by 



(Aii)c (iii2>G 
(x2X\)g {^2)% 



aG _ ( \^l/G \^^^2Hj \ 



where 



(^i)g = (^2)g = J J dxidX2x]p^p{xi,X2) 



= i((T + T2), (47) 



and 



{xiX2)g = {X2Xi)g 



j j dx\dx2X\X2p%p{x\,X2) 



From Eq. (43), (44), (47) and (48) we can observe that the condition 



(48) 



(49) 



must be satisfied for having A^^ = A^^. The Gaussian distribution which satisfies this condi- 
tion is denoted by p^p{xi ,X2). 

The differential entropy of this Gaussian distribution is given by 



HljyNFixuX2)] = log2 2ne 



0A2 



(50) 



The differential entropy for pfif{xi ,X2) reads 

H[pf^p{xuX2)] = -II dxidX2PNF{xi,X2)\og2{pNF{xi,X2)) 



logy(\/2naAiP) ; — — / / dxidx2 It 

V2naAiPJJ { 

o{xi -X2Y 



X sint' 



4p2 

2(il -X2f 



\Qg2e' 



-logjsint 



(51) 



If r = and s 



'2^ 



, we find that 



H\p^fjp{xx,x2)] = \og2{V2noAiP) 



21n2 



1 - — I ds sint^5^1n sint^i^ 
Ai 



(52) 



The integral 



I~ I ^fisint^i^ln sint^s^, 



must be calculated numerically. It gives that / ~ —0.692. Then, the differential entropy of the 
near-field joint distribution is 

//[p|,^(ii,f2)] «log2(V2^AiP) + 1.434. (53) 

Therefore, the negentropy of the near-field joint distribution is 

N\pnf(x\,X2)] = H[p'^p{xi,X2)]-H[plp{xi,X2)] 

« log2 {iTte^J^^P^ -\og2{V2TtaAxP) - 1.434 

« log2 ^\/27rey^^ - 1.434 

« 0.22. (54) 



Thus, the total non-Gaussianity of SPDC is 

nG^« 0.15 + 0.22 = 0.37. (55) 

D. Further properties of nG^ 

Property 1. iiG^ = iff the spatial two-photon state of SPDC is a Gaussian state. 
Proof. Since 

nG^ = Nnf +Nff = 0, (56) 
and that negentropy is always nonnegative, we have that N[pfip{xi,X2)] ~ N[ppp{qi,q2)] — 
0. Thus, for the joint probabilities we have that H[ppp{qi,q2)] = H[ppp{qi,q2)] and 
H[PNFi^uX2)] ^H[pf;p{xi,X2)], such that 

dqidq2PFF{qi,q2)iog2PFF{qi,q2) = 

J dqidq2Pff{qi,q2)iog2Ppp{quq2), (57) 

and 

dxidX2PNFi^l,X2)iog2PNF{^UX2) = 

dxidx2pffp{xi,X2)\og2PNF{xi,X2)- (58) 
Since the differential entropy // is a continuous and monotonic function, it holds that 

p9f (^1 > ^2)l0g2 Pff (^b = 

PFF{qiA2)^Og2PFF{q\A2), (59) 

and 

Pnf {^\ ,X2)\og2 Pnf {x\,X2) = 

PNF{x\,X2)\og2PNF{x\,X2). (60) 

Then, necessarily we have that 

PFFi^Um)^ PFF{^i,^2), (61) 

and 

PNpi^uh) PNFi^U^2)- (62) 

Thus, we conclude that the two-photon state is a Gaussian state. 

In the opposite direction we consider that the two-photon state is a Gaussian state. Then, the 
joint probability density functions of transverse position and momentum are Gaussian func- 
tions. So, from the definition of negentropy we must have that 

N[pipixi,X2)] =N[pfpiqi,q2)]^0. (63) 

So, the total non-Gaussianity of the spatial two-photon state becomes 

nG^ = 0. (64) 



Property 2. If the spatial two-photon state is a product state then nG is additive, i.e., nG = 
2 X nG^. 

Proof. For a set of two random and independent variables, the joint probability density func- 
tions are given by the product of the probability density functions associated to each variable, 
i.e., p{^i , 1^2) = )/'(&)■ Such type of probability density function has the covariance matrix 
Aij defined by the elements 

An = (^f)-(^i)' 

^1^2>-(^1>(&> 

J d^2^2p{^2)-j d^i^ipi^,) j d^2^2p{^2) 



=1 SI 



Pi^l) 



A22 



A12 = 



(65) 



(66) 



A21 = 



0, 
A12 



0, 



(67) 
(68) 



and a expected value vector jU = {(i^i),(^2)}- 

The Gaussian distribution p^{£,i,£,2) with the same covariance matrix and expected value 
vector is given by 



p''{^lA2) = 



1 



27r\/detA 



exp 



1 



(69) 



(70) 



\/27rAii \Jl%h-22 
where^ ={<^i,^2}- 

Note that these marginal Gaussian distributions have the same variance and expected value 
that and p{£,2). Since the random variables are independent, it holds that 

H[p{^U^2)]^H[pi^,)]+H[p{^2)], 
and consequently, we have that the negentropy of p{^i , ^2) can be written as 

N[p{^U^2)] = H[p^{^u^2)]-H[p{^u^2)] 

= //[/(<^i)]-^/[/^(^i)]+^b'^(&)]-^[p(&)] 

= N[pm+Nlp{M- 

If the two-mode spatial state of SPDC is a product state, the near- and far-field prob- 
ability density functions are written as Pff{qi,q2) ~ /'f fl^Off f(?2) and pfjp{xi,X2) = 
PNFi^i)PNFi^2)- Therefore, the properties mentioned above hold. Then, the total non- 
Gaussianity of the spatial two-photon state of SPDC will be given by 

"^^ = Nh4Fiqhq2)]+Nljyf,p{xuX2)] 

= N[Pff{9i)] +N[p'FFm] +N[p'^p{x,)] +N[pi,ix2)] 
= nG^'^+nG^^ (72) 



(71) 



nG' 



where we have used that nG^' = N\plp{qi)] + N[plp{xi)] and nG'^'- = A^[/4f (92)] + 

Due to the symmetry present in the two-photon wave function [see Eq. (3) and Eq. (4) of the 
main paper], we have that nG*^' = nG*'^ = nG^ . Then 

«G^ = 2 X nG*', (73) 

which proves our statement. 

Property 3. The non-Gaiissianity of the spatial state of SPDC decreases under partial trace, 
such that nG^ > nG^. 

Proof. To demonstrate this we first study the negentropies of the far- and near-field marginal 
distributions, while considering two limiting cases: (i) P ^ 1 and (ii) P very large. The joint 
distributions of the far- and near-field planes are given by Eq. (27) and Eq. (40), respectively. 
Thus, the marginal distributions are given by 



PFpih) =y dqiPFpiquqi), (74) 

and 

PNpixi) J dxipffp{xi,X2). (75) 
Doing the following change of variables: u = j{qi— qi) and v = ^^jT^' obtain that 

PfFfe)°= J due'7^'-"^'''^^^\mc^u^, (76) 

and 



V 2 



dve-^' V+'-f) siniV. ill) 



pirixi) °= J 

Considering the limiting case (i) we may write the marginal distributions as 

PFF{q2)°^ J du8{u+Pq2)sinc^u^ oc smc^{P^q2), (78) 

and 

p^f(x2)c<e-24. (79) 

In such case, it is clear that the negentropy of f^/rfe) is null [N[pfip(x2)] = 0]. The negentropy 
of the far- field marginal distribution will be the same of the function sine ^(P^ ^2)- It gives that 

NlpUiqi)] « log, - 0-4447 = 0.154. (80) 

Then, for values of P < 1, nG'^ w 0.154 such that nG^ > nG'^ . 
In the other limiting case (ii), we have that 

(81) 

and 



PNFixi)' 



11 ■ 1-^1 

smfv" o= sint -4r. 

p2 



(82) 



Now, the negentropy of p^ppiqi) is null {N\p^pp{qi)] = 0]. The negentropy of the near-field 



marginal distribution will be the same of the function sint^^. It gives that 



N\pIp{x2) 




(83) 



Thus, for larger values of P, nG^ w 0.224 which is also smaller than nG^ . 

Intermediate cases, for values of 0.002 <P <3, have been calculated numerically. The values 
of the negentropies of the far- and near-field marginal distributions are shown in Figure 3(b) of 
the main paper In Figure 3(d), iiG^ is plotted in terms of P. As one can see, in such intermediate 
cases, one always have that nG^ > nG^ . For P ~3 the value of nG^ w 0. 175. For values of P > 
9, nG^ already tends to its maximal value which is nG'^^,^ — 0.224 (See Figure 5 below). This 
proves that the non-Gaussianity of the spatial state of SPDC decreases under partial trace. 



0.23 
0.225 

0.22 
0.215 

Z 0.205 
0.2 
0.195 
0.19 



10 



Fig. 5. Here we show the non-Gaussianity of the marginal distributions for larger values of 
P. 



E. Calculating the non-Gaussianity of the state of the spatially correlated down- 
converted photons using the QRE 

For comparison purposes, we calculated the value of the non-Gaussianity measure based on the 
QRE, 5b [27,28]. This was done by assuming the case of perfect phase-matching to simplify 
the calculation of the mean values in the covariance matrix. Then, we considered that no corre- 
lations are present between the position of one photon and the momentum of the other photon. 
Therefore, the cross joint probabilities p{xi,qj) (where i,j = 1,2 and / 7^ j) are described by a 
product of the marginal position and momentum distributions p{xi) and p{qj)- In such case the 
covariance matrix can be written, considering the order V (i 1 , ^ 1 , f 2 , ^2 ) , as 



V 







11 







i 1 



73 



(84) 



The Gaussian state associated to this covariance matrix has a purity of /i = 0.44 [28], and its 
Von-Neumann entropy reads 5 = 1.08. Since the state of Eq (1) is pure, we have that 5b = 1.08, 
which does not depend on the value of P. 

F. nG^ at the limit wlien P < 1 

Here we discuss the behavior of nG^ when considering the limit of P <C 1, and that one of the 
down-converted photon is detected in transverse points around the origin. First we analyze the 
negentropy of far-field joint distribution: 



(85) 



In the limits considered, the differential entropies can be written as H[p{^i,^2)] ~ 
H[p{£,i\£,2)] + H[p{£,2)], such that the negentropy of the far-field joint distribution may also 
be rewritten as 



Nip^ppiquqi)] = 



H\p^p{qi\q2)]+H[pfp{q2)] 
-H[plp{q,\q2)]-H\i4p{q2)] 
H[p3Aqi\q2)]-H[p'ppiqi\q2)] 
+N[j7'pp{q2)]. 



(86) 



Note that Ppf{qi\q2) does not have the same expected value and variance of the conditional 
probability density function Pff{qi\q2)- Thus, in general, we can not define H[p^{qi\q2)] — 
H[ppp{qi 1^2)] as the negentropy of the conditional far-field distribution. 

Considering the joint far-field distribution (27), the conditional distribution when q2—Q will 
be given by 

plriqMi = 'i) °^ e-^ &mc^' 



4 



(87) 



When P <C 1, the sine function is much larger than the Gaussian function, such that we can 
approximate this distribution as 



PFFi.qi\q2 = Q) 



1 



2n 



Therefore, the differential entropy is 

^[PfF(^ife = 0)]«^log2(27re) 



(88) 



(89) 



Note that Eq. (88) is a normal distribution with the expected value equal to zero and the variance 
equal to 1 . 

Now let's consider the Gaussian distribution p'^p (^1 , ^2) of Eq. (3 1). It gives the conditional 
Gaussian distribution Ppp{qi \q2) at ^2 = equal to 



PFFiqi\q2 = o) = 



VJ+p^ 



-rr 



6n 



(90) 



In the case when P <C 1, this Gaussian distribution has also the expected value equal to 
zero and the variance equal to 1. Thus, it corresponds to the Gaussian distribution with the 
same expected value and variance of Ppp{qi\q2 = 0). Hence, in this limit, we can define the 
negentropy of the far-field conditional distribution in terms of the Ppp{qi \q2 = 0), such that 

N[p'ppiq,,q2)] ^N[l4p{qi\q2)] +N[,4piq2)]. (91) 

Now we consider the negentropy of the near-field joint distribution (40). If f ^ 1, then the 
near-field conditional distribution is given by 

pLi^i \X2 = 0) « ^ sint2 J^, (92) 

where we assumed (7=1 for simplicity. 
The variance of this function is 

of 9 <; 4A2P^ 

(Mi|i,=o) = J dxixip%p{xi\x2 = 0) = , (93) 

where A i and A2 are the constants defined above. Now we consider the near-field Gaussian joint 
distribution p^f{xi,X2)- The near-field conditional distribution associated with this Gaussian 
approximation, when P ^ 1 , is given by 



^-(^^l^^^0)^2;^^"^' '''' 

As one can note, p'^p (xi \x2 = 0) has the same expected value and variance of pfjp {x\ \x2 = 0). 
Therefore, in this limit we can write the negentropy of the near-field joint distribution as 

NIj4p{xuX2)] =NIj4p{xi\x2)] +N[p'np{x2)1 (95) 

where N[pf/p{xi \x2)] is defined in terms of the differential entropy of p'^p{xi \x2). 

Thus, when we consider P <C 1 we can write the total non-Gaussianity of the spatial two- 
photon state of SPDC as 

nG^w«G'^ + «G*' (iffP< 1), (96) 

where nG^ = N[pf,p{xi \x2)] +N[p^pp{qi \q2)], and nG'^ = N[pj,p{qj)] +N[pfjp{xj)], with ; = 
1,2. 
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